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Contrary to a widely accepted view the phase diagrams of La2-a:Sr3;Cu04 and YBa2Cu306+y, 
in spite of similarities are remarkably different. Both the electric conduction properties and the 
commensurate/incommensurate spin ordering properties differ dramatically. It is argued that the 
role of disorder in YBCO is insignificant while the bilayer structure is crucial. On the other hand 
in LSCO the intrinsic disorder to a large extent drives the properties of the system. The developed 
approach explains the low-temperature magnetic properties of the systems. The most important 
point is the difference with respect to the incommensurate spin ordering, including the difference 
in the incommensurate pitches. The present analysis demonstrates that the superconductivity is 
intimately related to the incommensurate spin ordering. 



In early days of high temperature superconductiv- 
ity there was a belief that the phase diagram of 
La2_a;Sr3.Cu04 (LSCO) represents a generic phase dia- 
gram of cuprate superconductors. Nowadays it has be- 
come clear that, in spite of similarities, there are very 
important differences between different cuprates. LSCO 
and YBa2Cu306+y (YBSO) are the best experimentally 
studied compounds in the low doping regime. This is why 
the present work addresses these compounds. In LSCO 
the doping level of Cu02-planes p practically coincides 
with Sr concentration, p x, while in YBCO, because of 
the partial filling of oxygen chains, the doping level is dif- 
ferent from the oxygen concentration y. In LSCO doping 
gives way to superconductivity at p > Psc ~ 0.055 and in 
YBCO at p > Psc ~ 0.065, see Fig. [H At first sight this 
indicates full similarity. However, I will argue that the 
mechanisms behind psc in those two compounds arc dif- 
ferent and the closeness of the two values of psc is purely 
accidental. An important observation is that the normal 
state electrical resistivities at p < p^c are very much dif- 
ferent. At low temperature, T < lOOA', and at doping be- 
low the superconductivity threshold the in-plane resistiv- 
ity of LSCO exhibits [J, Q the Mott variable-range hop- 
ping regime p cx expKTg/T)^/"^}. This indicates strong 
localization of holes in the Neel and the spin glass regions 
of the LSCO phase diagram. These arc the regions la and 
lb in FigUJ On the other hand, the in-plane resistivity in 
YBCO at p < Psc shows only logarithmic dependence on 
temperature, p oc ln(C/T), indicating weak- localization 
regime 0, 01 ■ This is the region 1 on the YBCO phase 
diagram, Fig[T] For example, at p sa 0.04 the in-plane 
resistivity of LSCO is about 5 times larger than that of 
YBCO at T = 10/i , and the same ratio is about 10'^ at 
T ~ IK. Thus, role of disorder in LSCO below the su- 
perconductivity threshold is crucial, while in YBCO the 
disorder is a relatively minor issue. 

The magnetic properties of the compounds are also 
very much different. The three-dimensional antiferro- 
magnetic (AF) Neel order in LSCO disappears at doping 
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FIG. 1: Schematic low-doping and low-temperature phase di- 
agrams of LSCO and YBCO. 

LSCO: la.AF order coexists with diagonal incommensurate 
spin structure, strong localization of holes, lb. Diagonal in- 
commensurate spin structure, strong localization of holes. 

2. Parallel incommensurate quasistatic spin structure, super- 
conductivity. 3. Parallel incommensurate dynamic spin struc- 
ture, superconductivity. 

YBCO: l.AF order, weak localization of holes. 2. Parallel in- 
commensurate quasistatic spin structure, superconductivity. 

3. Parallel incommensurate dynamic spin structure, supercon- 
ductivity. 



p « 0.02 and gives way to the so-called spin glass phase. 
The incommensurate magnetic order has been observed 
at low temperature in neutron scattering. This order 
manifests itself as a scattering peak shifted with respect 
to the AF position. The incommensurate scattering has 
been observed even in the Neel phase where it coexists 
with the commensurate one. In the Neel phase, the in- 
commensurability is almost doping-independent and di- 
rected along the orthorhombic b axis 0] . In the spin-glass 
phase, the shift is directed along the b axis, and scales 
linearly with doping In the underdoped supercon- 
ducting region (0.055 ^ p ^ 0.12), the shift still scales 
linearly with doping, but it is directed along one of the 
crystal axes of the tetragonal lattice In YBCO the 
commensurate three-dimensional AF order exists up to 
p w 0.065, see Fig. [TJ Moreover, there are indications 
that there is a narrow window around this doping where 
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superconductivity and the commensurate AF order co- 
exist 0. Recently the incommensurate quasistatic spin 
ordering along the tetragonal a* direction has been ob- 
served within the superconducting phase of YBCO (9l.[loj 
at doing p « 0.085. The ordering becomes fully dynamic 
above p « 0.1 [lH. Last, but not least, the observed in- 
commensurate wave vector in YBCO at p « 0.085 (ol. [loj 
is of a factor two smaller than the incommensurate wave 
vector in LSCO at the same doping Q . On the other 
hand, ai p ~ 0.12 the incommensurate wave vectors in 
LSCO and YBCO are equal 0,[il|- 

The phase diagram of underdoped LSCO has been ex- 



plained in Refs 
is, to a large extent 



13, 13, Physics of this compound 
driven by disorder. At low tem- 



perature each hole is trapped in a hydrogen-like bound 
state near the corresponding Sr ion. Each bound state 
creates a spiral distortion of the spin background. The 
distortion is observed in neutron scattering. So, the state 
at 0.02 < p < 0.055 is not a simple spin glass, it is a 
disordered spin spiral. Both the lower and the upper 
boundaries of this region are determined by the size of 
the bound state. The upper boundary, p = 0.055, is 
a percolation point of isolated bound states. After the 
percolation the superconductivity becomes possible, and 
simultaneously direction of the spin spiral must rotate by 
45°. The rotation is driven by the Pauli principle. The 
role of disorder at p > 0.055 is only marginal. Here the 
spin-spiral state suggested long time ago by Shraiman 
and Siggia [l3| is realized. Most importantly, the state 
is superconducting and the spin spiral becomes dynamic 
at p> 0.12 0. 

The present work is aimed at underdoped YBCO 
where, according to data on conductivity, the role of dis- 
order is practically insignificant. (This is consistent with 
the fact that a diagonal spin structure has been never ob- 
served in YBCO.) The following two issues are addressed. 
1) Why does the AF order survive up to a pretty large 
hole concentration p « 0.06 — 0.07? 2) Why is the pitch 
of the incommensurate spin order different from that in 
LSCO? It will be demonstrated that both these issues arc 
closely related and they are due to interlayer hopping. 

The present analysis of the single Cu02 -layer is based 
on the two-dimensional t — t' — t" — J model at small 
doping. After integrating out the high energy fluctua- 
tions one comes to the effective low energy action of the 
model [l^. Importantly, the integration of the high en- 
ergy fluctuations is a fully controlled procedure, the small 
parameter justifying the procedure is the doping level, 
p <C 1. The effective low-energy Lagrangian is written in 
terms of the bosonic n- field {n^ = 1) that describes the 
staggered component of the copper spins, and in terms 
of fermionic holons ip. I use the term "holon" instead 
of "hole" because spin and charge are to large extent 
separated, see The holon has a pseudospin that 

originates from two sublattices, so the fermionic field ip 
is a spinor acting on pseudospin. Minimums of the holon 



dispersion are at the nodal points qo = (±7r/2, ±7r/2). 
So, there are holons of two types (= two flavors) corre- 
sponding to two pockets. The dispersion in a pocket 
is somewhat anisotropic, but for simplicity let us use 
here the isotropic approximation, e (p) ^/^P^ , where 
p = q— qo. The lattice spacing is set to be equal to unity, 
3.81A 
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All in all, the effective Lagrangian reads 
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The first two terms in the Lagrangian represent the usual 
nonlinear a model. The magnetic susceptibility and the 
spin stiffness are x± ~ 0.53/8 w 0.066 and p^ w 0.18 [i3|. 
Hereafter the antiferromagnetic exchange of the initial t- 
J model is set to be equal to unity, J sa 130 meV 1. 
Note that ps is the bare spin stiffness, therefore by def- 
inition it is independent of doping. The rest of the La- 
grangian in Eq. ^ represents the fermionic holon field 
and its interaction with the n- field. The index a = 1,2 
(flavor) indicates the pocket in which the holon resides. 
The pseudospin operator is ^a, and = (1/v^, ±l/v^) 
is a unit vector orthogonal to the face of the MBZ where 
the holon is located. A very important point is that the 
argument of in Eq. ([1]) is a "long" (covariant) momen- 
tum, V = — i V + ^f? • [n'x Vn] . An even more important 
point is that the time derivatives that stay in the kinetic 
energy of the fermionic field are also "long" (covariant), 
J)t = dt + ■ [ft X n] . While the scmiclassical behaviour 
is determined by the Shraiman-Siggia term (the last term 
in ID)), the covariant derivatives arc crucial for quantum 
fluctuations and in particular for stability of the system. 

The effective Lagrangian ^ is valid regardless of 
whether the n-field is static or dynamic. In other words, 
it does not matter if the ground state expectation value of 
the staggered field is nonzero, (n) ^ 0, or zero, (n) = 0. 
The only condition for validity of ([T]) is that all dy- 
namic fluctuations of the n-field are sufficiently slow. 
The typical energy of the n-field dynamic fluctuations 
is E cross oc p^''^, see Ref. and it must be small com- 
pared to the holon Fermi energy ep (x p. The inequality 
Ecross ^ is valid up to optimal doping, p w 0.15. So, 
this is the regime where ([Ij is parametrically justified. 
Numerical calculations within the t — t' — t" — J model 
with physical values of hopping matrix elements give the 
following values of the coupling constant and the inverse 
mass, g w 1, /? Ri 2.2. On the other hand the fit of the 
neutron scattering data on LSCO gives g ~ 1, (3 ~ 2.7, 
which is in good agreement with the t — t' — t" — J model, 
see discussion in Ref. [lij . 

The dimensionless parameter 
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plays the defining role in the theory [l6|. If A < 1, the 
ground state corresponding to the Lagrangian ([T]) is the 
usual Neel state and it stays coUinear at any small dop- 
ing. If 1 < A < 2, the Neel state is unstable at arbitrarily 
small doping and the ground state is a static or dynamic 
spin spiral. Whether the spin spiral is static or dynamic 
depends on the doping level. The pitch of the spiral is 
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If A > 2, the system is unstable with respect to phase sep- 
aration and/or charge-density-wave formation and hence 
the effective long-wave-length Lagrangian ([T|) becomes 
meaningless. By the way, the pure t — J model {t' — 
t" = 0) is unstable since it corresponds to A > 2. Using 
values of g and /3 found from fit of experimental data, 
one obtains that for LSCO A ~ 1.30. 

How the described above physics is changed in case of 
YBCO? Due to the bilaycr structure the magnon spec- 
trum in YBCO is split into acoustic and optic mode 3] ■ 
The optical gap is about 70meV. This is substantially 
smaller than the maximum magnon energy ^ 2J ~ 
260meV. Therefore, the bilayer structure cannot sub- 
stantially influence values of the effective coupling con- 
stant g and the inverse mass (3 which are due to magnetic 
fluctuations with the typical energy scale ~ 2 J. So, one 
should expect that values of these parameters in YBCO 
are close to that in LSCO. The holon dispersion in YBCO 
is split into bilaycr bonding and antibonding branches 
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The splitting A is most likely due to the hole hopping via 
the interlayer oxygen chain sites. In any case both the 
LDA calculation [l§| and the ARPES measurements [2^ 
indicate the band splitting at nodal points about A ~ 
lOOmeV. In the present work A will be used as a fitting 
parameter. The splitting A brings additional nontrivial 
physics in the system. 

Let us impose the coplanar spiral configuration on the 
system ni — (cosq-r, sinq-r, 0) , 7?2 = — (cosq-r, sinq- 
r, 0) , where q is directed along the CuO bond [q cx (1,0) 
or q C)C (0, 1)]. Here ni and n2 correspond to the two lay- 
ers. Note, that fii and n2 remain antiparralell at any 
given point r, hence there is no an admixture of the op- 
tic magnon to the ground state configuration. The single 
holon energy spectrum is shown schematically in Fig. [2l 
There is the bonding- antibonding splitting A, and within 
each band there is a splitting between different pseu- 
dospin states, ±gg, because of the spiral. Populations 
of the four bands shown in Fig. [2] depend on the doping 
level p and on the spiral wave vector q. When calculating 
energy, one has to remember that there are two planes. 
Therefore, the elastic energy per unit area is 2 x psg^/2, 
and the hole density per unit area is 2p. It is convenient 




FIG. 2: Schematic dispersion of a holon. b,a -corresponds to 
bonding and antibonding brandies, and ± corresponds to the 
spiral splitting. 



to define the following characteristic concentration 
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A straighforward calculation shows that the dependences 
of the band filling and energy on doping p and pitch q are 
the following. Only the filled bands arc noted below. In 
the case p < po/2 the both fe-bands are filled at g < ^p 

and at q > —p only the band 6_ is filled. The energy is 
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In the case < P < Po the bands h± are filled at 

q < qi, the bands b± and a_ are filled at < q < q2, 
and b-, arc filled at g > 52- The energy reads 
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In the case p > po the bands b±, a± arc filled at q < q^, 
the bands b± and a_ are filled at 93 < 9 < 92, and 5_, 
a_ arc filled at g > (72- The energy reads 
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In these formulas 91 

and93 = f (f-f). 

The minimum of the energy with respect to 9 gives 
the equilibrium spiral pitch Q at a given doping level 
p. The result depends on A. For 1 < A < | the pitch 
stays zero for p < po, then for po < p < pi, where pi = 
0.5po/(A — 1), the pitch is 
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FIG. 3: Incommensurate pitch versus doping. 
Left: the regime 1 < A < 1.5. The solid line is the theory pre- 
diction for YBCO, the parameters are A = 1.3, po = 0.065. 
The dashed line is the theory prediction for LSCO. The red 
circles represent the YBCO neutron scattering data from 
Refs. [gl, [lol. [Till . The blue squares represent the LSCO neu- 
tron scattering data from Ref. 

Right: theoretical prediction for the pitch in the case 1.5 < 
A < 2. 



and finally at p > pi the pitch is given by the single layer 
formula ([3]). This behaviour is illustrated in Fig.[3l^Left). 
In the case | < A < 2 the spin spiral pitch stays zero un- 



til the critical concentration Pc = t 1 
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and then it jumps to the single layer value ([3]), see 
Fig. [JljRight). I would like to reiterate once more that 
the considered picture is valid for both the static and 
the dynamic spirals. Ultimately, the spiral always be- 
comes dynamic at p > 0.1 — 0.12, see Ref. [l^. Clearly 
at I < A < 2 the jump at p = Pc is the first order phase 
transition. On the other hand, po and pi at 1 < A < | 
are Lifshitz points. 

The value A w 1.3 has been obtained from the fit of 
LSCO data [l^. The parameter A cannot be influenced 
by the relatively weak interlayer coupling, therefore the 



same value should be used for YBCO [2l|. According 
to Refs. 0, Q the AF order in YBCO extends up to 
p = 0.06 - 0.07, so let us take po « 0.065. This is 
another parameter of the theory. Having these two pa- 
rameters one can predict the incommensurate wave vec- 
tor in YBCO. The prediction is shown in Figl3{Left) by 
the solid line The theory agrees very well with neutron 
scattering data shown by red circles 0, [l^, [lH. Using 
([5]) one finds the value of the bonding-antibonding split- 
ting, A w 70mcV. This is consistent with LDA calcu- 
lations [li] and with ARPES data [13 • In the same 
FiglHjLeft) the single layer theoretical Q(p) is shown 
by the dashed line, and the neutron scattering LSCO 
data 01 are shown by the blue squares. 

The developed theory is based on the small-p expan- 
sion. Therefore, it is not surprising that at p > 0.12 the 
experimental data start to deviate from the theory. Note 
also, that the single layer formula ^ is not applicable to 
LSCO at p < 0.055. The region p < 0.055 in LSCO corre- 
sponds to the strong localization regime and the relevant 
theory was developed in Ref. 



erty of underdopcd cuprates. However, the ordering 
properties and the phase diagrams of the single-layer 
LSCO and of the double-layer YBCO are remarkably 
different, see Fig. [TJ It is shown that while in LSCO 
the intrinsic disorder to a large extent drives the mag- 
netic properties, the role of disorder in YBCO is practi- 
cally insignificant while the bilayer structure is crucial. 
The present analysis demonstrates also that the super- 
conductivity is intimately related to the incommensurate 
spin ordering. In LSCO this relation is masked by the 
intrinsic disorder, superconductivity is impossible in the 
strongly-localized regime and therefore p^c is determined 
by percolation. However, in YBSO the correlation be- 
tween superconductivity and incommensurate spin order- 
ing is clear, the critical concentration for onset of super- 
conductivity practically coincides with that for onset of 
the incommensurate spin order. Physical mechanisms 
behind this observation will be considered elsewhere. 

I am grateful to O.K. Andersen and V. Hinkov for im- 
portant discussions. 
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